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Recent scanning tunneling microscopy measurements on cuprate superconductors liave revealed 
remarkable spatial inhomogeneities in the single-particle energy gap. Using cellular dynamical mean- 
field theory, we study the zero temperature superconducting properties of a single-band Hubbard 
model with a spatial modulation of the electron density. We find that the inhomogeneity in the 
electronic structure results in a substantial spatial variation in the superconducting order parameter 
and single-particle energy gap, reminiscent of the experimental results. In particular, we find that 
the order parameter and gap amplitudes in the hole-rich regions are significantly enhanced over the 
corresponding quantities in a uniform system, if the hole-rich regions are embedded in regions with 
smaller hole density. 
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I. INTRODUCTION 

The field of strongly correlated electron systems has 
been one of the most active areas of materials science 
for decadesii In particular, transition-metal oxides in- 
cluding high- Tc cuprates^ and colossal magnetoresistancc 
manganites are fascinating both from fundamental sci- 
ence and practical application points of view. Recent ex- 
perimental studies revealed that various transition-metal 
oxides exhibit spatial inhomogeneity in the electronic 
structurei^ For the cuprates, such inhomogeneity has 
been reported to have random, checkerboard or stripe- 
type spatial modulationi^"— 

Although its origin and physical properties remain a 
center of debate, it is clear that spatial inhomogene- 
ity in the cuprates greatly affects the superconducting 
properties^"— One of the most intriguing experimental 
observations is the positive correlation between the su- 
perconducting gap and the oxygen dopant atoms^i and 
the anticorrelation with the distance between a Cu atom 
and an apical O atom, d^, in the scanning tunneling mi- 
croscopy (STM) expcrimentsjii Based on an early mean- 
field study of an inhomogencous tJ model, Wang et alt^ 
conjectured that the inhomogeneity arises from varia- 
tions in the local oxygen dopant level. Other papers have 
addressed the inhomogeneity by introducing variations 
in coupling constants on a phenomenological mean-field 
level ji^ii^ More recently, Mori et ali^ proposed a link be- 
tween the variation in pairing interactions and the posi- 
tion of out-of-plane oxygen atoms, that involves the cova- 
lency between the in-plane dx2-y2 band and the apical p 
orbitals, and the screening of the in-plane Coulomb repul- 
sion on the Cu orbitals by the apical oxygen polarization. 
Here we study this problem in a microscopic model with 
inhomogeneity in the electronic structure using a non- 
perturbative approach that goes beyond a static mean- 
field treatment. 

In particular, we investigate the effect of spatial inho- 



mogeneity in the electronic structure on the supercon- 
ducting (SC) properties of a Hubbard model using cel- 
lular dynamical mean- field theory (CDMFT) Jlii^ Moti- 
vated by the result in Wang et alt^ and Mori et alt^ 
where out-of plane dopant atoms affect the in-plane hole 
density, we consider a two-dimensional Hubbard model 
with planar density modulations. We find that the d- 
wave gap amplitude is significantly increased in regions 
with increased hole density, if they are surrounded by re- 
gions with smaller hole density. The spatial modulation 
of the gap is found to become as large as a factor 2 and 
goes well beyond the variation one would expect from 
the difference in hole density. The SC order parameter 
also shows similar variation. The enhanced d-wave gap in 
the overdoped region is consistent with the experimental 
reports if the main effect of oxygen atoms is to locally 
increase the hole concentrationAii 



II. MODEL AND METHOD 

We consider a two-dimensional (2D) Hubbard model 
for the Cu d electrons: 

H ^^Ein^cr-t ^ {dl^dja + H.c.)+U^n,^n,i. (1) 



Here, di^ is the annihilation operator of an electron with 
spin a at site i, ni„ = d\^dia, t the nearest-neighbor 
transfer integral, and U the on-site repulsive Coulomb 
interaction. Here, only the nearest-neighbor hopping t is 
included. Ei is a site-dependent on-site potential, leading 
to a spatial inhomogeneity in the hole density. We con- 
sider two types of potential modulations as illustrated in 
Fig. [TJ A checkerboard-type variation in the local po- 
tential where Si alternates between Si = -\-V and —V on 
2x2 plaquettes and a stripe-like variation in which the 
modulation is unidirectional. 

We study the zero temperature SC properties of the 
uniform and inhomogencous systems using the CDMFT. 
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FIG. 1: Schematics of the potential modulation. Left: 
checkerboard-type modulation and right: stripe- type mod- 
ulation. On the sites belonging to a 2 x 2 plaquette denoted 
by -(-(—), the on-site potential is -\-V { — V). 



For a uniform system, the CDMFT maps the bulk lat- 
tice onto a cluster embedded in a dynamic mean-field 
that describes the rest of the system and is determined 
self-consistently. Here, we use a small 2x2 cluster. This 
approximation may be justified in the strong-coupling 
[U) limit, in which the coherence length becomes suf- 
ficiently small so that the superconducting behavior is 
well approximated on a 2x2 plaquette. In addition, 
the small cluster size allows us to solve the effective 
cluster problem using the Lanczos exact diagonalization 
(ED) techniqu o^^i^° to directly access dynamical quanti- 
ties such as the density of states without additional ana- 
lytic continuation. 

For the inhomogeneous system, we map the lattice onto 
a number Nc of in-equivalent 2x2 clusters, indicated by 
the bold lines in Fig. [1] where iVc is the number of clus- 
ters per unit cell. In our case we have Nc = 2. In this 
approximation, the self-energy depends on the cluster in- 
dex Ri and the intracluster coordinates fi. This means 
that dynamic correlations are retained within a cluster, 
but neglected between the clustersi ^^i^^ The clusters are 
then linked by the self-consistency condition. This ap- 
proach can be regarded as a cluster extension to the 
DMFT study of stripes in Ref. H^. The choice of clusters 
may not be unique. For example, one may take clusters 
comprised of two sites of the +V region and two sites 
of the —V region. However, in this case, the local sym- 
metry inside the cluster is severely broken. We expect 
this to give rise to a strong suppression of d-wave pairing 
inside the cluster. Hence this choice would not allow us 
to explore the physics we set out to explore, such as the 
modulation of a local pairing gap through the proximity 
to a region with different doping. 

With our choice of clusters, the impurity Hamiltonian 
is given by 

Himpj ^ Hi Eial^^aa^^ + ^jw^ctac,.,^ + /i.c.j 

ai(T aia 
+ Al^.yCr{aa2aaa3-a- + aaiattal-cr) + ft-.C.j. (2) 

Here, Hj represents the lattice Hamiltonian, Eq. ([T|), re- 
stricted to the Ith region (/ = l,...,Nc). We intro- 
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FIG. 2: (Color online) The superconducting order parameter 
'I' as a function of carrier density A'^ for the 2D Hubbard 
model computed by CDMFT with the ED impurity solver at 
T = 0. Black symbols are the results for the uniform systems. 
Colored (shaded) symbols connected by a thin line denote the 
results with a checkerboard-type potential modulation (upper 
panel) and a stripe- type potential modulation (lower panel). 



duced two kinds of bath orbitals a{= 1,2), which de- 
scribe the rest of the system and whose energy level and 
hybridization strength with a cluster site are given by 
Ea and Uq,, respectively. Thus, our cluster model con- 
sists of four interacting orbitals and eight bath orbitals. 
We also introduced two parameters representing the su- 
perconducting correlations for bath a, Aax and Aay. For 
the uniform system and the checkerboard-type potential, 
the local tetragonal symmetry is retained and therefore 
Aay = — Aax for the d-wave pairing. On the other hand 



in the stripe-type potential, the tetragonal symmetry is 
broken and Aax and Aay are independent parameters. 
Even in this case, the difference between Aay and — Aax 
was found to be very small. 

The self-consistency condition of the ED CDMFT is 
closed using the conjugate gradient minimization algo- 
rithm with a distance function defined on an imaginary 
frequency axis with the effective inverse temperature 
/3 = 50/t, and short cut-off frequency oJc = 1.5t^ We 
set the Coulomb repulsion to U = IQt and note that our 
results are insensitive to this choice as long as U is suffi- 
ciently large. Potential modulation V will be at most on 
the order of t, and Hubbard bands do not overlap with 
the low-energy quasiparticle part. Therefore, the effect 
of the short cut-off frequency is expected to be small. 
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TABLE I: Potential modulation V for the checkerboard-type 
pattern, the resulting carrier densities and A''_ , and super- 
conducting order parameters ^Z^. fj, is the chemical potential. 
Results with V = correspond to the uniform system. In the 
actual computations, we fix ^ -I- F and vary fj, — V. 
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TABLE II: Potential modulation V for the stripe-type pat- 
tern, the resulting carrier densities N+ and N- , and supercon- 
ducting order parameters ^'^^(y)- A* is the chemical potential. 
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III. RESULTS 

Numerical results for the d-wave SC order parameter 
"if — "i/x = ^'^y (^C ~ (citCi-i-54,)) as a function of carrier 
concentration A^ {J2a "-o-)ciuster for the uniform system 
are shown as black dots in Fig. [5] These results are in 
excellent agreement with the previous work by Kancharla 
et al^ We find an optimal carrier density A^°p* ^0.85 
which maximizes ^ and the superconductivity vanishes 
for densities smaller than A^ ^ 0.76. 

We now turn to the results for the inhomogencous sys- 
tems for the checkerboard and stripe-like modulations as 
shown in Fig. [TJ We plot the order parameter of clus- 
ter J, at the corresponding densities N = Nj = 
a ''^'^) cluster,! in Fig. [2] as colored (shaded) symbols 
connected by thin lines. One clearly observes that in 
the cluster with lower electron filling (larger hole density) 
is significantly enhanced over that in the uniform system 
at the same filling. For a range of fillings, it becomes even 
larger than the maximum value for the uniform system. 
At the same time, ^ in the cluster with the larger filling 
is reduced relative to the homogeneous system. Potential 
modulation, chemical potential, and the resulting carrier 
densities and SC order parameters are summarized in 
Table ID Here, / = -|-(— ) corresponds to the cluster in 
which the potential is increased (decreased) by V, thus 
< N^. In the actual computations, we fix fi + V and 
vary fjL — V . One can also see that V ^ {) recovers the 
uniform system. 

This enhancement of superconductivity is similar to 
that observed in our previous study of superlattices 
made up of underdoped and overdopcd Hubbard layersi^ 
There wc argued that the enhancement in the overdopcd 
layers was due to their proximity to the larger pairing 



scale in the underdoped layers. The similarity between 
the two results and the nature of approximation lead us 
to speculate that the origin of the increased order param- 
eter in the present study has similar origin. 

To support this conjecture, we plot the results for the 
system with the stripe-like modulation in the lower panel 
of Fig. [5] The corresponding numerical data are listed 
in Table [Hi Here, we plot both 'J'^: and —'^y How- 
ever, the difference between them is at most '^1% and 
therefore the corresponding points are almost indistin- 
guishable (see also Table HT]). As in the checkerboard 
case, the order parameter is enhanced in the cluster with 
larger hole doping. However, it is clear that the en- 
hancement is significantly weaker than in the former case 
with checkerboard-type potential. For the checkerboard 
modulation, the maximum 5* is about 30 % larger than 
the maximum value in the uniform system, while in the 
striped case the enhancement is about 10 %. We believe 
that this difference comes from the fact that, for the 
checkerboard modulation, the hole-rich cluster is fully 
surrounded by hole-poor clusters, which have a larger 
pairing scalei^ On the other hand, in the striped case, 
the hole-rich cluster is in proximity to only two hole-poor 
clusters and also two hole-rich clusters, which reduces the 
effect. In both cases, 4' in the hole-poor clusters is de- 
creased significantly by the coupling with the hole-rich 
clusters. The reduction is stronger in the stripe-type po- 
tential. This fact may partly be explained by the break- 
ing of the local tetragonal symmetry. These results sug- 
gest that the SC order parameter is quite sensitive to the 
local environment. 

It is also interesting to discuss the dependence of the 
order parameters in the hole-rich, ^f"*", and hole-poor re- 
gions, XP", with increasing potential or equivalently, 
decreasing filling A^-|_ and , respectively. With increas- 
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FIG. 3: (Color online) Local DOS of the Hubbard model 
with the checkerboard-type potential modulation. Each panel 
shows the local DOS in the clusters with smaller and larger 
for a given potential modulation. Down (up) triangles in- 
dicate the low-energy structure for the smaller (larger) A'' re- 
gion. To compute the self-energy, a small imaginary number 
77 — 0.125t was added to the real frequency lu. 



ing V and decreasing -/V+, follows the same trend as 
the order parameter in the homogeneous system, i.e., it 
first increases, then decreases with decreasing N+. It is 
enhanced over the order parameter in the uniform system 
through the proximity of the hole-rich region to the larger 
pairing scale in the hole-poor region. An optimal value 
of V that maximizes results from a delicate balance 
between the proximity of the larger pairing scale from the 
hole-poor regions, and the reduction in the pairing scale 
with increasing doping in its own hole-rich region. The 
order parameter in the hole-poor region, in contrast, 
decreases monotonically with decreasing filling N~ . This 
may be explained by a strong reduction in the pairing in- 
teraction in the hole-poor region by its coupling to the 
smaller pairing scale in the hole-rich regions. 

It would be interesting to determine whether the en- 
hancement of pairing we find in the hole-rich regions, 
translates to an increase in Tc- Our zero temperature 
and small cluster calculations, however, do not allow us 
to speculate on this question. Finite temperature calcu- 
lations on larger clusters similar to those carried out for 
a striped Hubbard model in Ref. [2^ will be necessary to 
gain more insight into this question. Work along these 
lines is currently underway. 

From the variation in the SC order parameter, one 
may expect that the d-wave gap in the density of states 
(DOS) is also modulated in the inhomogeneous sys- 
tems. To see this effect, we computed the local DOS 
for different strengths of the checkerboard-type potential 
modulations. We calculated the fc-space average of the 
momentum-dependent lattice Green's function to obtain 
the local Green's function, and the local DOS is given as 
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FIG. 4: (Color online) d-wave gap amplitude A as a function 
of the average carrier density for the 2D Hubbard model 
computed by the CDMFT with the ED impurity solver at 
T = 0. Black symbols are the results for the uniform systems. 
Colored (shaded) symbols connected by a thin line denote the 
results with a checkerboard-type potential modulation. 



its imaginary part. The lattice Green's function has the 
same form as that used to close the self-consistency con- 
dition for the CDMFT but with the self-energy as a func- 
tion of real frequency. When computing the self-energy, 
we introduced a small imaginary number rj on the order 
of O.lt to the real frequency. We confirmed that different 
choices of ?/ give similar results. Results presented below 
were obtained using rj = 0.125i. Because of the finite 
number of orbitals in the ED CDMFT and finite broad- 
ening 77, we cannot access the very low-energy behavior. 
Yet, we believe that the main structure of the local DOS 
is accessible. 

Results of the local DOS are shown in Fig. [31 Each 
panel corresponds to a different strength of potential 
modulation, and the resulting electron densities in hole- 
rich and hole-poor clusters are indicated. One can see 
that the low-energy structure (peak or shoulder) for the 
smaller region is located at higher energy than for the 
larger A^ region. 

In order to see the dependence of the gap amplitude A 
on the carrier density modulation more clearly, we plot 
A as a function of the average density A^ in Fig. |4] for 
the checkerboard-type potential modulations. Here, we 
defined A as the half distance between the lowest energy 
peaks or shoulders above and below a; = 0. For compari- 
son, we also plot A for the uniform system as black dots. 
The monotonic increase in A with approaching A^ = 1 is 
consistent with the previous reporti^ From this figure, 
one clearly sees that the gap amplitude in the hole-rich 
region is significantly enhanced over that in hole-poor re- 
gion, as well as that of the corresponding uniform system 
at the same hole density. We stress here that this result 
is robust against the choice of broadening parameter rj. 

From Figs. [3] and m we also notice that the effect of a 
potential modulation on the gap amplitude in the under- 
doped region is weaker than that on ^f; A is reduced less 
than when the potential modulation exists. We argue 
that this is caused by a near cancellation between the re- 
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duced order parameter ^' and the enhanced quasiparticle 
weight Z both due to the presence of the hole-rich region 
nearby as follows. The d-wave order parameter ^ is pro- 
portional to the anomalous self-energy at a; — > 0, and the 
real part of the self-energy is given by (1 — 1/Z)a;. The 
combination between the two effects leads to the gap am- 
plitude A oc . This insensitivity may become impor- 
tant to understand the doping-depcndent gap amplitude 
in the presence of inhomogeneity. 

We emphasize that the enhancement of the c?-wave gap 
in the hole-rich region relative to that in the hole-poor 
region is quite striking and opposite to the behavior in 
the uniform system. It is also contrary to the previous 
theoretical result based on a static mean-field theory.— 



IV. DISCUSSION AND SUMMARY 

Finally, let us briefly discuss the possible implication 
of the present study on the STM results. In Ref. H, it 
was reported that the gap amplitude and the position of 
dopant oxygen atoms have positive correlation. Further- 
more, m Ref.[nl an anticorrelation was reported between 
the gap magnitude and the distance between the Cu and 
apical O, dyi. If the main effect of the dopant oxygen and 
shorter (La is to increase the local hole concentration, as 
reported in the calculation by Mori et al.,— our obser- 
vation that the region with larger hole density (smaller 
filling N) has a larger gap is consistent with the experi- 
mental reports. 

Estimating potential modulation from experimental 
data is not easy because it requires precise position of 
charged (dopant) ions and the electron-density modula- 
tion. But, taking the experimental uncertainty of Ref@, 
density modulation is expected to be less than 10 %. 
In our calculations, the potential modulation V < t is 
enough to induce such a small density modulation. This 
V (< 0.5 eV) is smaller than the one simply estimated 
by assuming point dopant-charged ions. For example, 
in La2-KSra;Cu04, potential modulation due to the A- 
site dopant ions becomes of order of 1-2 eV with the 
in-plane lattice constant 4 A and the dielectric constant 
e ~ 3.5.— In reality, other effects such as the screening 
due to the other conduction electrons reduce the mod- 
ulation. Therefore, our potential modulation must be 
regarded as the self-consistently determined one by in- 
cluding screening effects, and V t seems to be in the 
realistic range. We would like to emphasize that our re- 
sults show that even a very weak density modulation, on 
the order of 5 % and therefore well within the experimen- 
tal uncertainty, is enough to result in a significant, i.e., 
factor 2, variation in the local gap amplitude. 

However, the variation in the gap magnitude we ob- 
serve of about a factor of two is less than the factor 3 
that has been reported in Refs. and [ll|. Certainly, 
a single-band Hubbard model is too simple to fully re- 
solve this issue, and additional degrees of freedom could 
play an important role, such as the apical oxygen as pro- 



posed in Ref. [l^ One possible scenario to resolve the 
present discrepancy is the screening effect by the oxy- 
gen ions. This reduces the on-site Coulomb interaction, 
resulting in a further enhancement of the superconduct- 
ing order parameter and the d-wave gap within the 
strong-coupling regimci^ 

As reported in the STM measurements in Ref. [^, both 
the average gap amplitude and the width and asymmetry 
of the distribution of gap magnitudes increase as the dop- 
ing decreases, with long tails on the large gap side of the 
distribution. Assuming that the observed inhomogeneity 
is related to a carrier-density modulation, our study can 
provide a possible explanation. Based on our results, we 
can speculate that the large gap tail in the underdoped 
samples originates from hole-rich islands with large gaps 
residing in the hole-poor environment. The strong asym- 
metry is naturally expected because the volume fraction 
of the large-gap hole-rich regions is small because of the 
small average hole concentration. 

Recently, Parker and coworkersi^ studied the gap- 
amplitude modulation in the over-doped cuprates 
Bi2Sr2CaCu208+5 as a function of temperature. It was 
revealed that medium-sized gap regions surrounded by 
large-gap regions persist to higher temperatures than re- 
gions with similar gap sizes in the proximity of small- 
gap regions. We speculate that both the large-gap and 
medium-gap regions are over-doped and the former is 
closer to the optimal doping. Thus, the proximity effect 
is weaker compared with that between under-doped re- 
gions and over-doped regions as we studied in this paper. 
To confirm this conjecture, it is necessary to make a link 
between the gap amplitude and the local carrier density. 
It is also desirable to perform the similar experimen- 
tal analysis using the under-doped and near-optimally- 
doped cuprates. 

To summarize, we have studied the superconducting 
properties of an inhomogeneous two-dimensional Hub- 
bard model for the high- Tc cuprates with modulations 
in the electronic hole density using a cellular dynamical 
mean-field approximation. We have found that the mod- 
ulation in hole density results in a considerable spatial 
variation in the superconducting order parameter which 
is quite sensitive to the local environment. In particu- 
lar, the order parameter in the hole-rich regions is sig- 
nificantly enhanced over the order parameter in a corre- 
sponding uniform system, if the hole-rich region is em- 
bedded in a region with smaller hole-density. Further- 
more, we have found that the d-wavc gap amplitude in 
the density of states shows similar modulation, with an 
enhanced gap in the hole-rich regions. These results are 
reminiscent of the variation in gap amplitudes reported in 
recent STM measurements of cuprate superconductors. 
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